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Abstract. In this paper we study an optimal control problem associated to 
a linear degenerate elliptic equation with mixed boundary conditions. The 
equations of this type can exhibit the Lavrentiefl phenomenon and non- 
uniqueness of weak solutions. We adopt the weight function as a control in 
L^{Q). Using the direct method in the Calculus of variations, we discuss the 
solvability of this optimal control problem in the class of weak admissible 
solutions. 
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1. Introduction 

The aim of this work is to study the existence of optimal solutions in coeffi- 
cients associated to a linear degenerate elliptic equation with mixed boundary 
condition. By control variable we mean a weight coefficient in the main part 
of the elliptic operator. The precise answer of existence or none-existence of 
an L^-optimal solutions is heavily depending on the class of admissible con- 
trols. Here are the main questions: what is the right setting of optimal control 
problem with L^-controls in coefficients, and what is the right class of admis- 
sible solutions to the above problem? Using the direct method in the Calculus 
of variations, we discuss the solvability of this optimal control problem in the 
so-called class of weak admissible solutions. 

In this paper we deal with an optimal control problem in coefficients for 
boundary value problem of the form 

{-div {p{x)Vy) + y = f in 
y = on Yd, 

/'(^)S=0 on Tat, 

where / E L^(J7) is a given function, the boundary of J7 is made of two disjoint 
parts do. = F/) U rjv, and p is a measurable non negative weight on a bounded 
open domain Q. in M^. 

Several physical phenomena related to equilibrium of continuous media are 
modeled by this elliptic problem. In order to be able to handle with media 
which possibly are somewhere "perfect" insulators or "perfect" conductors (see 
[9]) we allow the weight p to vanish somewhere or to be unbounded. 

Though numerous papers (see, for instance, [H 13 El [IHl EZl [El [E] and 
references therein) are devoted to variational and non variational approaches 
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to problems related to ([T]), only few papers deal with optimal control problems 
for degenerate partial differential equations (see for example [21 Ej ) . This can 
be explained by several reasons. Firstly, boundary value problem ([T|) for locally 
integrable weight function p may exhibit the Lavrentieff phenomenon, the non- 
uniqueness of weak solutions, as well as other surprising consequences. So, in 
general, the mapping p i— > y{p) can be multivalued. Besides, the characteristic 
feature of this problem is the fact that for every admissible control function 
p with properties prescribed above, the weak solutions of ([1]) belong to the 
corresponding weighted Sobolev space W^^P{Q.,pdx). In addition, even if the 
original elliptic equation is non-degenerate, i.e. admissible controls p are such 
that p{x) > a > 0, the majority of optimal control problems in coefficients have 
no solution (see for instance [llj). 

The optimal control problem we consider in this paper is closely related to 
the optimal reinforcement of an elastic membrane [3j. Reinforcing an elastic 
structure subjected to a given load is a problem which arises in several appli- 
cations. The literature on the topic is very wide; for a clear description of the 
problem from a mechanical point of view and a related bibliography we refer 
for instance to the beautiful paper by Villaggio [14j . 

In the simplest case when we have an elastic membrane occupying a domain $7 
and subjected to a given exterior load / G L^(0), the shape u of the membrane 
in the equilibrium configuration is characterized as the solution of the partial 
differential equation 

— div [p{x)Vy^ + y = f in Q 

together with the corresponding Dirichlet and Neumann boundary conditions 
on 50. The reinforcement of the membrane is usually performed by the addition 
of suitable stiffeners, whose total amount is prescribed. Mathematically, this is 
described by a nonnegative coefficient p(x) which acts in and is associated 
with some weight coefficient in the main part of elliptic operator. As a result, the 
problem of finding an optimal reinforcement for the membrane then consists in 
the determination of a weight p{x) > which optimizes a given cost functional. 
In contrast to the paper [3] , we do not restrict of our analysis to the particular 
case of the reinforcement problems. We also do not make use of any relaxations 
for the original optimal control problem. 

2. Notation and Preliminaries 

Let be a bounded open subset of (A'^ > 1) with a Lipschitz boundary. 
We assume that the boundary of is made of two disjoint parts 

with Dirichlet boundary conditions on T/j, and Neumann boundary conditions 
on Ftv. Let xe be the characteristic function of a subset S C O, i.e. xe{x) = 1 
if x G i?, and XEi-c) = if x ^ E. 

LetC7^(M^;ri3) = C^(M^) : = on Td}. The space ^^'^O; Td) 
is the closure of C^{W'^ ■,r£)) in the classical Sobolev space W^'^{Q). For any 
subset E G U we denote by \E\ its iV-dimensional Lebesgue measure L'^{E). 
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Let p : — )■ M be a locally integrable function on such that p{x) > 
for a. e. X G M.^ . Then p gives rise to a measure on the measurable subsets 
of M'^ through integration. This measure will also be denoted by p. Thus 
p{E) = pdx for measurable sets E C . 

We will use the standard notation L^{Q,, pdx) for the set of measurable func- 
tions f on Q such that 

2, - 



We say that a locally integrable function p : — )• is a weight on Q if 

(2) p + p'^eLlM"")- 

Note that in this case the functions in L'^{Q,, pdx) are Lebesgue integrable on 

n. 

To each weight function p we may associate two weighted Sobolev spaces: 
Wp = W{n,pdx) and Hp = H{n,pdx), 
where Wp is the set of functions y G W^'^{Q; F/j) for which the norm 

(3) \\y\\p=(^Jjy^ + p\Vy\^)dxy^ 

is finite, and Hp is the closure of C^{Q;Td) in Wp. Note that due to the 
estimates 

(4) / \y\dx<( [ |ypdx)'^Vr/' < C\\y\\p, 



(5) J \Vy\dx<(^j \Vy\^pdxy^^(^J p'^ dxY^^ < C\\y\ 

the space Wp is complete with respect to the norm || • ||p. It is clear that 
Hp C Wp, and Wp, Hp are Hilbert spaces. If p is bounded between two positive 
constants, then it is easy to verify that Wp = Hp. However, for a "typical" 
weight p the space of smooth functions Cq°{CI) is not dense in Wp. Hence the 
identity Wp = Hp is not always valid (for the corresponding examples we refer 

to [a IIS]). 

Weak Compactness Criterion in L^{Q,). Throughout the paper we will often 
use the concepts of the weak and strong convergence in L^(il). Let {oejg^o 
be a sequence in L^{Q,). We recall that {ogj^^Q is called equi-integrable if for 
any 6 > there is r = t{6) such that Jg\ae\dx < 5 for every measurable 
subset S C of Lebesgue measure \S\ < r. Then the following assertions are 
equivalent: 

(i) a sequence {agj^^Q is weakly compact in L^(il); 

(ii) the sequence {agj^^Q is equi-integrable; 

(iii) given 6 > there exists A = X{6) such that sup / [agl dx < 5. 

s>0 J{\a^\>X} 

Theorem 2.1 (Lebesgue's Theorem). // a sequence {oej^^o ^ L^{^) is equi- 
integrable and Oe ^ a almost everywhere in then ^ a in L^{Q,). 
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Radon measures. By a nonnegative Radon measure on we mean a nonneg- 
ative Borel measure which is finite on every compact subset of Vl. The space 
of all nonnegative Radon measures on Q. will be denoted by M+(0). According 
to the Riesz theory, each Radon measure ^ G M+(i7) can be interpreted as 
element of the dual of the space Cq{^) of all continuous functions vanishing 
at infinity. Let M(yt; M^) denotes the space of all M^-valued Borel measures. 
Then// = (//i, . . . , ^uat) gM(J7;M^) <^ G C'q{VL), i = l,...,N. 

If is a nonnegative Radon measure on fi, we will use U{VL, dfi), 1 < r < oo, 
to denote the usual Lebesgue space with respect to the measure fj, with the 

corresponding norm ||/||L'(Q,d/x) = iln\fi^)\''dt^)^^''- 

Functions with Bounded Variation. Let / : ^ M be a function of L^{Q). 
Define 



where div (p = Yli=i 'BF'- According to the Radon-Nikodym theorem, if \Df\ < 
+00 then the distribution Df is a measure and there exist a vector-valued 
function V/ € [-L^(r2)]^ and a measure Dsf, singular with respect to the A^- 
dimensional Lebesgue measure restricted to Q, such that 



Definition 2.2. A function f £ L^{Q,) is said to have a bounded variation in 
ft if \Df\ < +00. By BV{Q) we denote the space of all functions in L^{Q) 
with bounded variation. 

Under the norm 



BV{0,) is a Banach space. It is well-known the following compactness result 
for i?y-functions: 

Proposition 2.3. The uniformly bounded sets in BV-norm are relatively com- 
pact in L^{Q,). 

Definition 2.4. A sequence {fk}'kLi C BV{Q) weakly converges to some f G 
BV{il), and we write fk ^ f iff the two following conditions hold: fk ^ f 
strongly in L^in), and Dfj, Df weakly* in M(Q;M^). 

In the proposition below we give a compactness result related to this conver- 
gence, together with the lower semicontinuity property (see [6]): 

Proposition 2.5. Let {fk}'kLi ^ sequence in BV{Q?j strongly converging to 
some f in L}{Q) and satisfying sup^^^j^ \Dfk\ < +00. Then 




\Df\ = sup 




Df = VfL''[Q + DJ. 




(i) : / G BVm and \Df\ < limmfk^^J^ \Dfk 

(ii) : fk^f in BV{n). 
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Convergence in variable spaces. Let {nk}keN^ be Radon measures such 
that /Ufc ^ /i in M+(r2), i.e., 

(6) hm / ipdnk= / ifdfi G Co(M^), 

where Co(M^) is the space of all compactly supported continuous functions. 
The typical example of such measures is 

dfJ'k = Pk{x) dx, dfi = p{x) dx, where < pk ^ p in L^{Q). 

Let us recall the definition and main properties of convergence in the variable 
L^-space (see [IS]). 

1. A sequence E d// A:)} is called bounded if lim sup / \vk\'^ dpk < 

A;— >oo Jq 

+ 00. 

2. A bounded sequence {vk G dpk)} converges weakly to v € L'^{Q,, dp) 



if 



lim / Vkfdpk = / vipdp V99 € C^{Q), 

'2 



and it is written as Vk ^ v in L (il, dpk). 
3. The strong convergence Vk ^ v m Lp'iyi, dpk) means that v € L^(il, dp) 
and 



(7) lim / VkZkdpk= I vzdp as Zk ^ z in L {Q.,dpk)- 

k-^"^ Jn Jn 

The following convergence properties in variable spaces hold: 

(a) Compactness: if a sequence is bounded in L^(r2, dpk), then this sequence 
is compact in the sense of the weak convergence; 

(b) Lower semicontinuity: if Vk ^ v in L'^{Q,, dpk), then 

(8) liminf / \vk\'^ dpk > / dp; 

Jn Jn 

(c) Strong convergence: Vk ^ v ii and only if Vk ^ v in L'^{Q, dpk) and 

2 J,, _ / „,2 



(9) lim / \vk\ dpk = V dp. 

fe-i-oojn Jn 

3. Setting of the Optimal Control Problem 

Let m € M+ be some positive value, and let ^1, ^2 be given elements of L^{Q) 
satisfying the conditions 

(10) ^1(2;) < 6(^) a.e. in n, ^ G L^fl). 

To introduce the class of admissible i?y-controls we adopt the following concept: 

Definition 3.1. We say that a nonnegative weight p is an admissible control 
to the boundary value problem 

(11) -div {p{x)Vy) + y = f m Q, 

dy 

(12) y = OonrD, p{x)—=OonTN, 
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(it is written as p ^ "^ad) 

(13) p^BV{Q), / pdx = m, ^i{x) < p{x) < £,2{x) a.e. in Q. 

Jn 

Here f G L'^{Q,) is a given function. 

Hereinafter we assume that the set Had is nonempty. 



Remark 3.2. In view of the property ()T0|) . we have the boundary value problem 
for the degenerate elliptic equation. It means that for some admissible controls 
p G Jlad the boundary value problem pip - ()12p can exhibit the Lavrentieff 
phenomenon, the nonuniqueness of the weak solutions as well as other surprising 
consequences. 

The optimal control problem we consider in this paper is to minimize the dis- 
crepancy between a given distribution yd G L^(ri) and the solution of boundary 
valued problem (fTT]) - (fT2]) by choosing an appropriate weight function p £ Jiad- 
More precisely, we are concerned with the following optimal control problem 
(14) 

I JQ Jn Jo. ) 



subject to the constraints (jTT|l - ()T3]l . 

Definition 3.3. We say that a function y = y{p, f) G Wp is a weak solution to 
the boundary value problem (jlip -(jl2 p for a fixed control p G Jlad if the integral 
identity 



fif dx 



(15) / [pVy ■V<f + y^)dx= / , 
Jn Jn 

holds for any ip G C^{Q;Te,)- 

It is clear that the question of uniqueness of a weak solution leads us to 
the problem of density of the subspace of smooth functions C^{Q; F/)) in Wp. 
However, as was indicated in [16], for a "typical" weight function p the subspace 
C^{Q] Tf)) is not dense in Wp, and hence there is no uniqueness of weak solu- 
tions (for more details and other types of solutions we refer to [H [lOl [151 [16]). 
Thus the mapping p y{p, f) is multivalued, in general. Taking this fact into 
account, we introduce the set 

(16) Ew = {{p, y) \ P & Had, y G Wp, (/>, y) are related by ([E]) } . 

Note that the set E]y is always nonempty. Indeed, let Vp be some intermediate 
space with Hp C Vp C Wp. We say that a function y = y{p, f) G is a 
V^-solution or variational solution to the boundary value problem (Ilip - (ll2p if 
the integral identity (I15p holds for every test function 99 G Ip. Hence, in this 
case the energy equality 

(17) / (|Vy|V + y') dx= [ fydx 

Jn Jn 

must be valid. Since 

fydx 



/, 

Jn 



< ( / fdxV^'l [ y'dxY^'<C\\y\ 



n ' ^Jn 
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for every fixed / € L'^{^1), it follows that the existence and uniqueness of a 
V^-solution are the direct consequence of the Riesz Representation Theorem. 



Thus every variational solution is also a weak solution to the problem (fTT 
(jl2p . Hence 'Ew ^ and therefore the corresponding minimization problem 



(18) inf_ I{p,y) 

is regular. In view of this, we adopt the following concept. 

Definition 3.4. We say that a pair {p^,y^) G L^(0) x W^'^{^;Ti:)) is a weak 
optimal solution to the problem (|13p -(|14 p if (p^,y^) is a minimizer for (jlSp . 

4. Existence of Weak Optimal Solutions 

Our prime interest in this section deals with the solvability of optimal control 
problem (|13p ~(jl4 p in the class of weak solutions. To begin with, we make 
use of the following results. Let {{pk,yk) € Sv^lfegN sequence of weak 

admissible solutions. 

Lemma 4.1. Let {pk}keN a sequence in 5iad such that p^ ^ p in L^(Q,). 
Then 

{Pk)~^ in the variable space L?{Q., pkdx). 



Proof. To prove this result we use some ideas of |16) . By the properties of the 
set of admissible controls "Jlad-, we have p'^^ < for every /c G N, hence the 
sequence {p^^j^jgpj is equi-integrable on il. Note that, up to a subsequence, we 
have Pk ^ p a.e. in VL. Since .^2^^ < p^^ < , Lebesgue Theorem implies 

(19) pl^^p-^\nL\n). 

Let if € C^{Vt) be a fixed function. Then the equality 

Pk^fPkdx= / Lpdx= I p^^ifpdx VA; € N 



n 



leads us to the weak convergence p^^ p^^ in L'^{i}, pkdx). Taking into 
account the strong convergence /?^^ in L^(il) and the fact that is a 

bounded domain, we get 

lim / \pk\^'^Pkdx= lim / p'^^ dx = / p^^ dx = / \p\~'^pdx. 

Hence, by the strong convergence criterium in the variable space L'^{d, pkdx), 
we come to the required conclusion. □ 

Our next step deals with the study of topological properties of the set of 
weak admissible solutions E]y to the problem pT]) - (fHp . 

Definition 4.2. A sequence {{pk,yk) ^ ^vfI^gn is called bounded if 



sup 

km 



\Pk\\BV{n) + WUkWl^in) + \\'^yk\\L^{Q,p^dx)N 



< +00. 



Definition 4.3. We say that a bounded sequence {{pk,yk) ^ ^v^lfceN of the 
weak admissible solutions t -converges to a pair {p,y) G BV{Q) x W^'^{Q) if 



8 GIUSEPPE BUTTAZZO PETER I. KOGUT 

(a) p in BV{Q); 

(d) Uk ^ y weakly in L?'{Q); 

(e) Vyfc ^ 9 L?{U,, pdx)^ in the variable space L'^{Q., pkdx)^ . 

Note that due to assumptions (fTO]l . (fT3l) . and estimates hke dl])-®, the 
inclusion y E VF^'^(il) is obvious. 

Lemma 4.4. Let {{pk,yk) £ ^lyjfcGN a bounded sequence. Then there is a 
pair{p,y) G BV{ft)xW^'^{^) such that, up to a subsequence, {pk,yk) (p^y) 
and y G Wp. 

Proof. By Proposition 12.31 and the compactness criterium of the weak conver- 
gence in variable spaces, there exist a subsequence of {{pk,yk) S ^vfI/cgNi 
still denoted by the same indices, and functions p G BV{Q), y G L'^{0,), and 
V G L'^{Q,pdx)^ such that 

(20) pk^pmL\n), 

(21) yk ^ y in L^iVl), ^yk ^ in the variable space L'^{ft,pk dx). 

Let us show that y G W^'^{0,), and v = Vy. Since ii < Pk < ^2 for every A; G N, 
([20|) and Lemma WA\ implv the property (see (fT9]) ) 

(22) pl^ ^ p-^ m L^{Vl), ^1 < p < ^2 a.e. in 

This yields that the sequence {VyfcjfcgN is weakly compact in L^{Q)^ . Indeed, 
the property of its equi-integrability immediately follows from the inequality 

/ \Vyk\dx<(f Pk^dxY^^if {Vykl^PkdxV^^ <C( f p^^ dxV^^ 

J E J E J ^ J E 

As a result, using the strong convergence {pk)~^ P~^ in the variable space 
L'^{Q, Pkdx) (see Lemma HTTI) and its properties, we obtain 

hm / Vyk -ipdx = lim / p].^Vyk ■ ipPk dx 



p ^v-ippdx= I V ■ il) dx 



'n Jn 

for all V G C^{n)^. Thus Vy^ ^ t; in L^{n)^. This implies that y G W^''^{n) 
and Vy = v. The inclusion y G Wp immediately follows from ()20p - (j2ip . □ 

Theorem 4.5. For every f G Ll^{R^) the set Ew is sequentially closed with 
respect to the t- convergence. 

Proof. Let {{pk,yk)}keN C be a bounded r-convergent sequence of weak 
admissible pairs to the optimal control problem (fTT]) ~ (fn|) . Let {po,yo) be its 
r-limit. Our aim is to prove that (To,yo) ^ '^w- By Lemma 14.41 we have 

(23) Pk ^ PO in L\n), po G BV{n), 6 < Po < 6 a.e. in 0. 

Then passing to the limit as /c — ?> oo in the relation pk dx = m, we conclude 
that Po S 3?ad) i-e. the limit weight function po is an admissible control. 

It remains to show that the pair (po,yo) is related by the integral identity 
(fT5]) for all ip G C^{^l]T£)). For every fixed A; G N we denote by {pk,yk) £ 
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BVioc{R^) X wl^^{M^) an extension of the functions {pk,yk) ^ '^w to the 
whole of space such that the sequence {{pk,yk)}keN satisfies the properties: 

(24) pk G BV{Q), ii<Pk< 6 a.e. in Q, 

(25) sup \\h\\BV{Q) + \\yk\\L^{Q) + \\^yk\\L^{Q,p^dx)t^ < +00 

ken L J 

for any bounded domain Q in M^. Hence, as done in Lemma [4.4l it can be proved 
that for every bounded domain Q C there exist functions po ^ BV{Q) and 
yo £ W^po such that 

(26) pk^PQ'va.L^{Q), ijk ^ Vq m. L^{Q), 

(27) V^A: ^ Vyo =5 L'^{Q,po dx)'^ in the variable space L?'{Q.,pkdx)^ . 
It is important to note that in this case we have 

(28) = yo and po = Po a.e. in 

In what follows, we rewrite the integral identity p5p in the equivalent form 
(29) 

{Vyk-V^Pk + ykV)XQ{x)dx= [ f^xnix)dx e C^{R^ ,Td), 



and pass to the limit in ([29]) as /c — )• oo. Using the properties (|26|) ~ (j27j) . and 
the fact that xn Xn strongly in the variable space L^{Q,pk dx), i.e. 

XnPk dx= xnPk dx — > l xnPo dx = XnPo dx 

JR'^ JR'^ Jr'^ 

we obtain 

I {S/yo-Vippo + yo^)xn{x)dx= [ fipxn{x)dx Mip G C^iR"" ,Td) 
Jr^ jrJv 

which is equivalent to 

/ (Vyo ■ Vippo + yoip) dx = ftpdx Vcp G C^(ri,r^). 
Jn Jn 

As a result, taking into account and the fact that yo € Wp^ (by Lemma 
14. 4|) . we conclude: yo is a weak solution to the boundary valued problem (flT]) - 
(fT^ under p = po- Thus the r-limit pair (ro,yo) belongs to set Hvi/, and this 
concludes the proof. □ 

We are now in a position to state the existence of weak optimal pairs to the 
problem pT]) - p^ . 



Theorem 4.6. Let € Lj^^iR'^) and ^,2 G Lj^^iR^) he such that < 6 a.e.in 
R^ and ^ G L]^^{R^). Let f G Lf^^(R^) and ya G L'^{n) be given functions, 
and assume that 3lad ^- Then the optimal control problem (fTT]) - (fH|) admits 
at least one weak solution 



opt . 
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Proof. Since the set of admissible controls Jlad is nonempty the minimization 
problem (fT8]l is regular (i.e. 7^ 0)- Let {{pk,yk) G ^w^lfceN be a minimizing 
sequence to (jlSp . From the inequality 



inf I{p, y) = lim 



\yk{x) - yd{x)\^ dx 



(30) 

+ / \^yk{x)\''pkdx + I \Dpk\ 
Jn Jn 

there is a constant C > such that 



< +00, 



snpWykh^n) <C, sup ||Vyfc 11^2 m ^^)iv < C, / \Dpk 
fcGN kern ^ ' Jn 



< C. 



Hence, in view of the definition of the class of admissible controls 3iad) the 
sequence {{pk,yk) G ^lyjfcGN is bounded in the sense of Definition 14.21 Hence, 
by Lemma [13] there exist functions p* G BV{Q,) and y* £ Wp* such that, up 
to a subsequence, {pk,yk) {p*^y*)- Since the set Hiy is sequentially closed 
with respect to the r-convergence (see Theorem 14. 5p . it follows that the r-limit 
pair (p*,y*) is an admissible weak solution to the optimal control problem (|lip - 
(i.e. {p*,y*) G Hiy). To conclude the proof it is enough to observe that 
by property ([8]) and Proposition 12.51 the cost functional / is sequentially lower 
r-semicontinuous. Thus 

I{p*,y*) <liramfl{pk,yk) = inf Iip,y). 

Hence {p*,y*) is a weak optimal pair, and we come to the required conclusion. 

□ 
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